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SUMMARY 

The generalized method of cells micromechanics model is utilized to analyze the tensile stress-strain response 
of a representative titanium matrix composite with weak interfacial bonding. The fiber/matrix interface is modeled 
through application of a displacement discontinuity between the fiber and matrix once a critical debonding stress has 
been exceeded. Unidirectional composites with loading parallel and perpendicular to the fibers are examined, as well 
as a cross-ply laminate. For each of the laminates studied, analytically obtained results are compared to experimental 
data. The application of residual stresses through a cool-down process was found to have a significant effect on the 
tensile response. For the unidirectional laminate with loading applied perpendicular to the fibers, fiber packing and 
fiber shape were shown to have a significant effect on the predicted tensile response. Furthermore, the interface was 
characterized through the use of semi-empirical parameters including an interfacial compliance and a “debond 
stress;” defined as the stress level across the interface which activates fiber/matrix debonding. The results in this 
paper demonstrate that if architectural factors are correctly accounted for and the interface is appropriately charac- 
terized, the macro-level composite behavior can be correctly predicted without modifying any of the fiber or matrix 
constituent properties. 


INTRODUCTION 

The usefulness of titanium matrix composites (TMCs) in the aircraft industry has been intensely studied. Conse- 
quently, the need for advanced analytical tools that predict the behavior of these materials has also increased. Par- 
ticularly important is the ability to model the effects of variations in the local microstructure and material 
composition on the overall, macroscopic response. 

Several factors have been found to play a significant role in the inelastic response (e.g., tensile and creep 
response) of TMCs; i.e., residual stresses, debonding and architectural effects to name a few. For example, Lerch 
and Saltzman (ref. 1) conducted an extensive series of experimental and analytical studies on the tensile deformation 
response of a silicon fiber reinforced TMC, that is, SiC/Ti-15-3. Specifically, when the material was loaded trans- 
verse to the fiber direction, the stiffness of the stress-strain curve decreased significantly once a certain stress level 
was reached. The decrease in stiffness was determined to be the result of interfacial debonding. Two factors were 
found to contribute to the debonding of the composite. First of ail, due to a cool-down from the heat treatment tem- 
perature, residual stresses were present in the fiber and matrix, which served to hold the constituents in place. These 
residual stresses needed to be overcome by mechanical loading in order for debonding to take place. However, 

Lerch and Saltzman found that the stress levels required for debonding were higher than what could be accounted 
for by residual stresses alone, indicating that there was also chemical bond strength holding the fiber and the matrix 
together. Therefore, both the residual stresses and chemical bond strength needed to be overcome by mechanical 
loading in order for debonding to take place. Expanding upon this study, Lerch et al. (ref. 2) conducted detailed 
finite element analysis of a [0/90] s SiC/Ti-15-3 TMC system, using gap elements to represent the fiber/matrix inter- 
face and including a cool-down in order to apply residual stresses. From the results of these analyses, the authors 
determined that interfacial debonding played a significant role in the tensile deformation of the composite, and that 
correctly incorporating the residual stresses was key in correctly predicting the onset of debonding. 

Majumdar and Newaz (ref. 3) also carried out an experimental investigation of a model TMC (SiC/Ti-15-3). 
Their studies found that when the composite was loaded perpendicular to the fiber direction, the tensile stress-strain 
curve had three distinct regions. First, an initial linear region occurred (denoted as Stage I, see fig. 7 in ref. 3), 
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during which the macroscopic behavior was elastic and the fiber and matrix remained well bonded. Second, once 
interfacial debonding occurs, a second linear region, with a lower stiffness than the original linear region (Stage IT) 
was observed. This decrease in stiffness was a result of fiber/matrix debonding; in that, once debonding took place, 
the stiff fibers had little contribution to the overall stiffness of the material. Finally, the macroscopic tensile response 
became nonlinear due to significant inelasticity in the matrix. The region was denoted as Stage III. Supplementary 
analyses also showed that for the composite examined in this work the debonding was primarily due to the overcom- 
ing of residual stresses by mechanical loading, with negligible bond strength. Further discussions on the effects of 
weak interfacial bonding on the transverse (perpendicular to the fiber) stiffness and strength of titanium matrix com- 
posites, including the effects of residual stresses can be found in reference 4. 

Robertson and Mall (ref. 5) conducted analytical studies using finite elements and a mechanics of materials 
based micromechanics approach on a TMC composed of Si C fibers in a TIMET AL 21 S 1 matrix. In the course of 
their analyses, the authors computed tensile stress-strain curves for loading perpendicular to the fiber direction and 
compared the results to experiments. Residual stresses were applied through a cool down process, and a weak fiber/ 
matrix interface was incorporated into their modeling. Upon examination of the results, the authors determined that 
the fiber/matrix bond was not due to residual stresses alone. Instead, they hypothesized that the fiber/matrix bonds 
possessed a finite strength, and debonding would not take place until the chemical or mechanical bonds failed due to 
the applied load. 

Arnold et aL, (ref. 6) examined in a qualitative sense the effects of fiber architecture and fiber/matrix bonding 
on the tensile and creep response of a TMC system (SiC/TIMETAL 2 IS). The analyses were carried out using the 
Generalized Method of Cells (CMC) micromechanics method (ref 7) and its framing computer code MAC (Micro- 
mechanics Analysis Code) (ref. 8). The analyses were carried out at elevated temperature; consequently the effects 
of residual stresses were neglected. The weak fiber/matrix interface was modeled as a distinct constituent with an 
elastic-perfectly viscoplastic constitutive law. The analyses indicated that both fiber architecture (i.e., the placement 
and spacing of the fibers within the composite) and interfacial bonding had a significant effect on the tensile and 
creep response of the composite, particularly when loaded perpendicular to the fiber direction. 

The current work 2 expands on the previous work of Arnold et aL, (ref. 6) in that it examines quantitatively the 
predictive capabilities of MAC (ref. 8) through the inclusion of residual stress effects and the comparison of experi- 
mental results. In the previous work isothermal analyses were conducted at one temperature level (650 °C) and 
residual stresses due to the cool-down history during heat treatment (processing) were not accounted for in the 
analysis. Here, the effects of residual stresses are shown to be quite significant with respect to capturing the onset of 
nonlinearity in the stress/strain curve (the “knee” in the curve) as well as the ultimate longitudinal strain to failure 
of the composite system. Furthermore, unlike the analyses conducted in (ref. 6), in the current work tensile stress/ 
strain response curves for [0], [90] and [0/90] laminates are compared to actual experimental results obtained at both 
room and elevated temperature. In this manner, a quantitative as well as qualitative measure of the effects of influ- 
encing factors (parameters) such as architecture, residual stresses and interfacial behavior on the tensile response 
can be examined. Additionally, the quantitative comparison of analytical to experimental results will demonstrate 
that by properly taking into account factors such as architecture, residual stresses and interface specifications, good 
correlation’s and predictions can be achieved with experimental data without modifying the fiber or matrix 
constitutive properties. 


ANALYSIS METHODOLOGY 
GMC 

The generalized method of cells (ref. 7), which is an extension of the micromechanics method known in the 
literature as the method of cells (ref. 10), is the micromechanics analysis method adopted throughout this study. In 
the original method of cells, a continuously (or discontinuously) reinforced composite is modeled as a doubly (or 
triply) periodic array of fibers or inclusions embedded in a matrix phase. The periodic character of the assemblage 
allows one to identify a repeating unit cell that can be used as a building block to construct the entire composite. The 
properties of the unit cell are thus representative of the properties of the entire assemblage. The unit cell consists of a 


! TIMETAL 21 S is a registered trademark of TIMET, Titanium Metals Corporation, Toronto, OH. 

2 This work was conducted back in 1996-1997 and is being presented at this time in order to provide expanded documentation of (ref. 9), 
and background/motivation for a future report (ref. 1 1). 
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single fiber subcell surrounded by three matrix subcells for continuous composites, and seven matrix subcells for 
discontinuous composites, hence the name method of cells. The Cartesian geometry of the repeating unit cell allows 
one to obtain an approximate solution for the stresses and strains in the individual subcells given some macroscopi- 
cally homogeneous state of stress or strain applied to the composite. The approximate solution to the thus posed 
boundary value problem is, in turn, used to determine the macroscopic (average) or effective properties of the com- 
posite and the effective stress-strain response in the inelastic region. 

In the generalized method of cells for continuous (or discontinuous) multiphased composites, the repeating unit 
cell is subdivided into an arbitrary number of subcells of phases. Specifically, for a two-dimensional analysis the 
unit cell is broken up into Np by Ny subcells, while for a three-dimensional analysis the unit cell is broken up into N a 
by No by Ny subcells (fig. 1). Therefore, this generalization extends the modeling capability of the original method 
of cells to include the following: (1) inelastic thermomechanical response of multiphased composites, (2) modeling 
of various fiber (phase) architectures including both shape and packing arrangements, (3) modeling of porosities and 
damage, and (4) the modeling of interfacial regions around inclusions, including interfacial degradation. 

The basic homogenization approach taken in the micromechanical analysis consists essentially of four steps. 
First, the representative volume element, RVE, (or repeating unit cell) of the periodic composite is identified. Sec- 
ond, the macroscopic or average stress and strain state in terms of the individual microscopic (subcell) stress and 
strain states is defined. Third, the continuity of tractions and displacements are imposed at the boundaries between 
constituents. For the case of imperfect fiber/matrix bonding, as will be described shortly a displacement discontinu- 
ity is applied between fiber and matrix subcells, which is related to the tractions across the interface. The above 
three steps, in conjunction with microequilibrium, establish the relationship between micro (subcell) total, thermal 
and inelastic strains and macro (composite) strains via the relevant concentration tensors. In the fourth and final step, 
the overall macro constitutive equations of the composite are determined. These four steps form the basis of the 
micro-to-macro mechanics analysis, which describes the behavior of heterogeneous media (refs. 12 and 13). The 
resulting micromechanical analysis establishes the overall (macro) behavior of the multiphased composite and is 
expressed as a constitutive relation between the average stress, strain, thermal strain and inelastic strain, in conjunc- 
tion with the effective elastic stiffness tensor. 

The resulting analytic constitutive law has been readily applied to investigate the behavior of various types of 
composites, given knowledge of the behavior of the individual phases. One advantage to this type of methodology is 
that any type of simple or combined loading (multiaxial states of stress) can be applied irrespective of whether or not 
symmetry exists. Similarly, one does not have to resort to different boundary condition application strategies. 
Another advantage of having available an analytical expression representing the macro elastic-thermoinelastic con- 
stitutive law, is that it ensures a reduction in computational costs and memory requirements when implementing this 
methodology, either in a stand-alone situation or within a finite element code (ref. 14). Evidence of this computa- 
tional savings (in excess of a thousand times) can be found in references 6 and 15, where a comparison between 
macroscopic stress-strain curves computed using GMC and FEA is made. There it was shown that for the same 
global response prediction, GMC was over three orders of magnitude faster in speed than the finite element method, 
and this was for simpler case of perfect bonding between fiber and matrix. In the case of weak interfacial bonding 
the difference in computational speed would increase even more with little or no loss in accuracy. 

To model imperfect fiber/matrix interfaces in metal matrix composites, two approaches are commonly utilized. 
One approach, which is used for example in reference 6, is to model the interface as a distinct constituent with its 
own material model. Another approach is to model the effects of interfacial debonding by applying displacement 
discontinuities between the fiber and matrix (ref. 16). In the case of finite element analyses, the latter approach is 
accomplished by applying gap elements between the fiber and matrix (ref. 2); whereas in MAC a modified 
Achenbach and Zhu (ref. 17) interface model has been developed and implemented and will be described in the next 
section. Furthermore the GMC formulation also admits a wide variety of physically based deformation and life con- 
stitutive models, with an example viscoplastic deformation model being discussed in the subsequent section. 


FTBER/MATRIX DEBONDING MODEL 

Weak interfacial bonding, due to the multilayer coating applied to the fiber prior to consolidation, is known to 
be present in TMCs (ref. 4) and significantly affects the transverse stiffness and strength of these composites. The 
fiber/matrix interface has been modeled in MAC by treating the interface as a distinct constituent (ref. 6). Modeling 
the interface in this manner allowed a consistent multiaxial representation during both loading and unloading of the 
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interface and the examination of the influence of initial imperfections and induced interfacial damage. However, 
for the current work, attempts to employ this method to simulate the interface resulted in exceeding long execution 
times for the analysis. Furthermore, the properties of this interface were unknown and needed to be arbitrarily 
assigned. 

In an attempt to improve the computational efficiency while still capturing the essence of the interfacial behav- 
ior, a new “elastic-perfectly plastic” flexible interface model, in the spirit of Achenbach and Zhu (ref. 17) and Jones 
and Whittier (ref. 1 8), was devised to model the interfacial behavior. In this method, a jump in the displacement 
field is applied between the fiber and matrix subcells when a critical debond stress is exceeded, while continuity of 
the traction vector is maintained. This critical debond stress activation concept was adopted since in at least some 
composites (ref. 5) the interface does not debond as soon as the stress across the interface becomes tensile (i.e., the 
mechanical clamping stress due to residual stresses is exceeded). Instead, due to chemical bond strength or other 
factors, the interface will only debond when a certain nonzero positive stress (“debond stress”) is applied across the 
interface. 

Such a concept is fully substantiated, by examining the elevated temperature tensile response of a TMC (SiC/ 
TIMET AL 2 IS) loaded perpendicular to the fiber direction (see the curve labeled G DB = 0 in fig. 7), where it is 
obvious that both the initial modulus and the knee in the stress-strain curve are significantly under predicted. These 
results indicated that debonding did not take place as soon as the stresses became tensile, implying that finite bond 
strength must have existed between the fiber and matrix for this composite. The interface model is described 
mathematically as follows: 


•/_[ 0 if of, < Gpg„ 

u n ~~ i _ . / . r / . 

[/?„ ct-„ if o n >a DBn 

/ f 0 if of < o DBt 
' 1 R r tj if oj>c Dm 


( 1 ) 


where 

u n displacement component normal to interface I 
u f displacement component traverse to interface I 
<5 n stress component normal to interface I 
<5 t shear stress component at interface I 
R n normal interfacial compliance 
R t tangential (shear) compliance 
a DBn normal debond strength 
®DBt tan g en ^ a ^ (shear) debond strength 

dotted quantities are time derivatives 

Note that the interfacial compliance (R n ) measured the degree of displacement discontinuity in the interface, which 
simulated the degree of debonding which was present. An R n of zero implied a perfect bond (no debonding), 
whereas an R n of infinity implied a fully debonded interface. R n was also defined by Aboudi (ref. 10) as the ratio of 
the interfacial thickness to the “Young’s modulus” of the interface. Thus, a R fJ of zero implied an infinite interfacial 
modulus, while a R ft of infinity implied an interfacial modulus of zero. The definition of and behavior of R r w as the 
same as for R with the only difference being in the physical interpretation of R t as the ratio of the thickness of the 
interface to its “shear modulus.” An important point to note is that in actuality the interface is a small area with 
undefined properties. Therefore, all of the material constants described here as being associated with the interface 
were in actuality semi-empirical parameters describing the physical mechanisms of debonding which were taking 
place. The exact processes used to determine the values of the parameters will be described in a later section of this 
report. Furthermore, in GMC the displacement continuity conditions could be applied to individual subcell bound- 
aries in specific coordinate directions. As will be discussed in greater detail in reference 11, this feature allowed for 
the simulation of partial interfacial debonding for a unit cell in which multiple subcells were used to model the fiber. 
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CONSTITUENT CONSTITUTIVE MODELS 


For this study, the material analyzed is a representative titanium matrix composite composed of continuous 
silicon carbide (SCS-6 3 ) fibers embedded in the titanium alloy matrix, TIMETAL 2 IS. The high-strength, high- 
stiffness, SCS-6 fibers were assumed in this study to be isotropic and linear elastic. Furthermore, as the variation of 
elastic properties with temperature was relatively small, the following material properties were utilized throughout, 
and were assumed to be independent of temperature: Young’s modulus, E = 400 GPa, Poisson’s ratio, v = 0.32; and 
coefficient of thermal expansion, a= 3.78X -6 °C~K 

TIMETAL 2 IS is a metastable beta titanium alloy, containing -21 percent alloying additions that has high 
strength as well as good creep and oxidation resistance. TIMETAL 21 S was developed for use in advanced metal 
matrix composites and its (isotropic) viscoplastic response has been characterized for the model of Bodner and 
coworkers (refs. 19 and 20), by Kroupa and Neu (ref. 21) as well as a generalized viscoplasticity with potential 
structure (GVTPS) model (ref. 22) by Arnold et al., (refs. 23 and 24). The GVIPS model of Arnold et al., (refs. 23 
and 24) was selected for use in the current study. In this model, specific forms of both the Gibb's and complemen- 
tary dissipation potentials were chosen such that a complete (i.e., fully associative) potential based multiaxial, 
nonisothermal, unified viscoplastic model was obtained. This model possesses a tensorial internal state variable and 
an evolutionary law that has nonlinear kinematic hardening and both thermal and strain induced recovery mecha- 
nisms. A unique aspect of this model is the inclusion of nonlinear hardening through the use of a compliance opera- 
tor, derived from the Gibb's potential, in the evolution law for the back stress. This nonlinear tensorial operator is 
significant in that it allows both the flow' and evolutionary laws to be fully associative (and therefore easily inte- 
grated), it greatly influences the multiaxial response under nonproportional load paths, and in the case of 
nonisothermal histories, it introduces an instantaneous thermal softening mechanism proportional to the rate of 
change in temperature. In addition to this nonlinear compliance operator, a new, consistent, potential preserving, 
internal strain unloading criterion was introduced to prevent abnormalities in the predicted stress-strain curves, dur- 
ing unloading and reversed loading of the external variables. These features make GVIPS accurate in relaxation, as 
well as able to handle any nonproportional loading that may be present in the composite. 

The flow and evolution equations were defined as follow's: 


0 





if F<0 
if F>0 


A 9 m 


Qijtnn^mnkl b kl 


b iJ 


if a^j < 0 
if Qijlij > 0 


a ij = Lijkl{Akl ~ ®klT)> 


( 2 ) 


where the dot represents time differentiation. A- is the deviatoric back strain tensor, a f j is the deviatoric back stress 
tensor, and T is temperature. Additional quantities were defined as follows: 


||ejy | = z\j t)j = (equivalent inelastic strain rate) 


^ijkl Qijkl 


bu = e •/ - 



/ 

3®i (p _ ^)G P ~ 2 aija kl 

3B 0 (\ + B lP GP- 1 ) 

K 2 0 [l + B l pG p - l (6p-5)\ j 

3pict^tf v m 3R a B 0 G q 1 

2k 2 0 4g k 2 0 J 

a ij 


(stiffness operator) 


( 3 ) 


0.. = -^-(l -+* B x pG p l )— ^(dynamic thermal recovery operator) 

J dT ' ' 


, 3tf/r 

K 


\SCS-6 is a registered trademark of Textron Systems Division. 
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where 


F = 

r = 

G = 

'i= 


'■in 


- Y ) (threshold function) 


/ 


- (yi 


(yield stress function) 


-4- (back stress function) 


3 3 

— , J 2 = — (stress invariants) 


( 4 ) 


and are elastic stiffness coefficients, H v [*] is the Heaviside unit step function, and <•> are Macauley brackets. 
The temperature-independent material parameters are: k q , n , Sj,/? and q, while the temperature -dependent material 
parameters are: K, ji o , Z? 0 , and Interpolation functions defined by Arnold et al., (ref. 23) were employed to deter- 
mine the material parameters for TIMET AL 21 S at temperatures other than the reference temperature of 650 °C. A 
limitation of the material characterization was that above 704 °C, material parameters had to be taken as those at 
704 °C. Material parameters at various temperatures are taken from (ref. 23) and presented in table I. 


COMPARISON OF ANALYTICAL COMPUTATIONS TO EXPERIMENTAL RESULTS 

The computationally efficient and comprehensive micromechanics analysis code, MAC, was utilized to obtain 
all of the results presented in this study. A full description of the usage and capabilities of MAC can be found in 
reference 8. Unless otherwise specified, a fiber volume fraction of 0.33, which was determined experimentally to be 
the nominal fiber volume fraction of the actual composite (ref. 23), and a pseudo-square idealization of the fiber 
shape 4 was utilized. A micrograph showing a portion of the actual specimen that was tested is shown in figure 2. As 
can be seen in the figure, while the distance between fiber layers (l.e., in the 2-direction, see fig. 3) is relatively con- 
stant, the distance between fiber centers within a layer (i.e., in the 3-direction, see fig. 3) varies significantly. The 
average ratio of the horizontal (3- axis) and vertical (2-axis) distances between fiber centers (R ratio) over the entire 
specimen was determined to be 1.11. Although, local fiber architecture can have a significant impact on the trans- 
verse material behavior (i.e., when a composite is loaded transverse to the fiber direction) it has little or no effect on 
the longitudinal response (i.e., when a composite is loaded in the fiber direction). Similarly, when a tensile load is 
applied perpendicular to the fiber direction, fiber/matrix interfacial debonding becomes significant and must be 
accounted for, whereas in the case of longitudinal loading (parallel with the fiber direction) fiber/matrix interface 
assumptions have no impact. The “elastic-perfectly plastic” flexible interface model discussed above was utilized to 
model the interfacial behavior, and will be discussed below. 

The behavior of [0], [90] and [0/90] laminates are examined in this report. Applying a longitudinal load parallel 
with the fiber axis (along the 1-axis in fig. 3) simulated the [0] laminate, whereas applying the load perpendicular 
(transverse, i.e., along the 2- or 3 -axes) to the fiber direction of a unidirectional laminate simulates a [90] laminate. 
For the [0] and [90] laminates, the doubly periodic, rectangular fiber array, GMC model was utilized (ref. 8). For the 
[0/90] laminate, the triple periodic, continuous reinforcement, GMC model was used as illustrated in figure 1 . 

As discussed earlier, residual stresses due to heat treatment have been found to significantly influence the onset 
of the nonlinear behavior of metal matrix composites. Consequently, in this study the residual stresses incurred dur- 
ing cool-down from heat-treatment (as all prior residuals due to consolidation are assumed to be annealed out) are 
considered. To incorporate the residual stresses, a uniform temperature cool-down profile is applied from the heat- 
treatment temperature to the temperature of interest. The starting temperature and cool-down rates were determined 
from the actual treatment profile. The results were examined at two temperatures of interest, i.e., 23 °C (room tem- 
perature, where residual stresses are significant) and 650 °C (elevated temperature, where residual stresses are 
negligible). In this way, the effects of residual stresses can be explicitly determined. Subsequent to the application of 

4 AIthough, the actual shape of the fiber is square, due to the averaging nature of the GMC formulation, the stress concentrations that one 
might expect at the comers are not present, therefore we use the term pseudo-square to remind the reader. 
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the cool-down profile, a monotonic strain-controlled, tensile test was simulated up to a strain level of 1.4 percent [0] 
or 2.0 percent for [90] laminates. The total applied strain rate for the [0] and [0/90] laminates was l.OxlO^/sec, 
whereas for the [90] laminate it was 1.667xl0^/sec. These total strain rates correspond to those utilized in the actual 
experiments (refs. 25 and 26). 


LONGITUDINAL TENSILE RESPONSE 

The simulated and experimental longitudinal tensile stress-strain results at elevated temperature are shown in 
figure 4, and results obtained at room temperature are shown in figure 5. In both figures, predictions made with and 
without the incorporation of residual stresses are compared to experimental results obtained by Castelli (ref. 27). 
Furthermore, in both figures the yield stress (c y ), defined as deviation from proportionality, is indicated in order to 
point out graphically the approximate point at which the curves become nonlinear. At elevated temperature, the in- 
clusion of residual stresses had a negligible effect (as expected) on the computed results, clearly indicating that the 
residual stress level was low at this elevated temperature, since it is near the heat treatment temperature. At room 
temperature, on the other hand, the residual stresses had a significant impact on the ‘yield stress” or the onset of 
nonlinearity (the “knee” in the curve) in the predicted stress/strain response. In figure 5, C R represents qualitatively 
the yield stress when residual stresses were incorporated into the simulation, and C^ R represents qualitatively the 
yield stress when residual stresses were not incorporated into the simulation. Incorporating the residual stresses into 
the simulation lowered the yield stress in the room temperature tensile stress-strain curve and provided a more favor- 
able match with the experimental results. These results indicated that at room temperature residual stresses were 
significant and needed to be incorporated into the simulation in order to make accurate predictions. 

For both temperature levels, at higher strain levels the analytical simulations predicted stresses slightly higher 
than observed in the experiments. While some of this discrepancy could be due to inaccuracies in the matrix consti- 
tutive model, it may also be due to fiber breakage occurring at higher strain (stress) levels in the experimental tests, 
which would lead to a softening of the experimental stress-strain curve. The fact that fiber breakage can take place 
before total failure of the material occurs is support for example by the work of Neu and Roman (ref. 28) and 
Brindley and Draper (ref. 29). 

In order to examine the issue of fiber breakage in more detail and given the fact that the longitudinal composite 
ultimate strength is dictated by the ultimate tensile strength of the SCS-6 fiber reinforcement, constituent fiber stress 
versus the total composite strain are superimposed on the overall stress-strain response of figure 5, (see fig. 6). 
Again, results are shown for cases with and without the incorporation of residual stresses. As the statistically signifi- 
cant ultimate tensile strength of the SCS-6 fiber reinforcement (ref. 30) varies between 2400 and 5170 MPa with the 
average being -4140 MPa, we have indicated in figure 6 potential fiber breakage levels (i.e., points A, B and C) 
starting at 3450 MPa and increasing at 350 MPa increments. Clearly the overall strain level at which these fiber 
stress levels occur is relatively close to the strain level at which the overall composite fails, with 3622 MPa being an 
accurate calibration estimate of the effective fiber strength, for both room and elevated temperature (not shown) 
calculations. Also an important observation relative to the inclusion of residual stresses due to processing can be 
gleaned from figure 6 in that when residual stresses are not applied, the composite strain level at which potential 
fiber breakage can occur (i.e.. A’, B' and C’) is significantly lower than the corresponding case with residual 
stresses. These results indicate that if residual stresses are not applied, the failure stress and strain of the overall 
composite could be significantly under predicted in an analysis. Consequently, given that the ultimate failure of the 
composite for this laminate and loading configuration is most likely due to fiber failure, the above results suggest 
that if an appropriate fiber breakage model was developed for MAC (see ref. 11) the analysis should accurately 
predict the ultimate longitudinal tensile failure strain of t he composite. Additionally, since the actual fiber tensile 
strength possesses a statistical distribution, see the vendor’s histogram (ref. 30), with its average at 4140 MPa, one 
would expect that if a statistical failure distribution was added to the fiber breakage model in such a way as to cap- 
ture the stress redistribution occurring due to the accumulation of fiber breaks a softer response (better representing 
the actual experimental results) would be predicted at high stress (strain) levels. Support for this idea can be seen in 
figure 6 wherein fiber breakage could initiate at composite stress levels as low as 2400 MPa at room temperature. 
The 2400 MPa fiber stress level corresponds to a composite strain level of 0.8 percent at which the experimental 
results become more nonlinear and start to deviate significantly from the computed results. This trend was observed 
both at room and elevated temperature. The results for both the overall composite and the fiber confirm that proper 
application of residual stresses can have a significant effect on accurately predicting both the material’s deformation 
behavior and ultimate tensile strength. 
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TRANSVERSE TENSILE RESPONSE 


The transverse (loading perpendicular to the fiber direction) tensile response for the [90] laminate was also 
simulated using MAC. To apply the transverse loads, strain controlled loading was applied along the 3-direction axis 
to a level of 2.0 percent at a strain rate of 1.667xlCT 4 /sec. The strain rate was chosen in order to facilitate compari- 
sons to available experimental data found in reference 26. The use of experimental data to correlate and validate the 
analytical response is the major difference between the current work and the previous work of Arnold et al., (ref. 6). 
In reference 6, the effects of architecture and the fiber/matrix interface bond were examined strictly from a qualita- 
tive viewpoint. Alternatively, in this work, the effects of the various parameters effecting transverse behavior (bond 
strength and architecture) are examined quantitatively through comparisons with experimental results. 


Correlation of Interfacial Parameters 

As discussed earlier, TMCs have been found to have a weak fiber/matrix interface and the effects of the inter- 
face become significant when tensile loads are applied perpendicular to the fiber direction. For the elastic-perfectly 
plastic interface model utilized, the debond strength, o DB , and the interfacial “compliance” need to be determined 
empirically. For this study, the composite was assumed to be perfectly bonded in the shear direction. Therefore, the 
tangential “compliance” R t was set to zero, and the tangential debond strength was set to a very high value. 

The first parameter to be correlated was the normal debond stress, o DS . For these analyses, the normal “com- 
pliance” R n was set to a relatively high value to simulate a fully debonded interface. Simulation results were exam- 
ined at both elevated temperature (650 °C) (fig. 7) and room temperature (fig. 8) (23 °C) and compared to exper- 
imental values (ref. 26). However, the elevated temperature results were the ones used for correlation 5 whereas, the 
room temperature results were used for verification. For both temperatures, the debond stress was varied from 
0 MPa (simulating zero bond strength) to 96.53 MPa. The results at elevated temperature (fig. 7) showed that when 
the debond stress was set to zero, the initial deviation from linearity (indicating the beginning of debonding) 
occurred almost immediately after application of the transverse tensile load, and at a much lower stress level than 
indicated by the experimental results. This was expected since minimal residual stresses were present at this tem- 
perature so that with zero bond strength, debonding would occur as soon as residual stresses were overcome. These 
results clearly indicate, consistent with reference 5, that finite bond strength exists for the SCS-6/TIMETAL 21 S 
composite independent from residual stresses. As the bond strength was increased, the stress level at which devia- 
tion from linearity began (again indicating the onset of debonding) increased. A debond stress of 96.53 MPa pro- 
vided the best correlation between the predicted and experimentally observed deviation from linearity. An 
interesting point to note is that varying the debond stress only appears to affect the stress level at which debonding 
commences while the remainder of the stress- strain response curves remain roughly parallel (similar). Considering 
the room temperature results, shown in figure 8, the initial deviation from linearity is clearly seen to be predicted 
reasonably well when the value of the debond stress is 96.53 MPa. Alternatively, when the debond stress was set to 
zero, the stress level at which debonding occurred (indicated by the onset of nonlinearity) was significantly under 
predicted. However, in sharp contrast to the elevated temperature response (fig. 7) debonding did not occur immedi- 
ately upon application of the tensile load, due to the higher level of compressive, thermally induced residual stress 
present at room temperature. These results appeared to confirm that for the SCS-6/TTMETAL 2 IS composite, both 
residual stresses and the finite bond strength needed to be overcome in order for debonding to occur at room 
temperature. 

Once the debond stress was correlated, the value of the normal interfacial compliance, R n was determined. Ana- 
lytical results were computed for both room and elevated temperature with the value of the interfacial compliance 
varying from 0 MPa' 1 (simulating a perfect bond) to a value of 4x10^ MPa -1 . Increasing the value of R n higher 
than 4X10" 4 MPa -1 did not change the results, which indicated that an R n value of 4X10 -4 MPa -1 simulated a fully 
debonded interface. The simulated and experimental (ref. 26) results are shown in figure 9 (elevated temperature) 
and figure 10 (room temperature). The results in both figures clearly demonstrate that when a perfect bond is 
assumed (R n equal to zero) the correlation between the experimental and analytical results is significantly over pre- 
dicted. Alternatively, simulating a fully debonded interface provided the best correlation to the experimental results. 


5 From longitudinal results, the residual stresses at elevated temeprature were determined to be negligible, consequently the initial deviation 
from linearity in the macro stress-strain can be attributed solely to the breaking of the fiber/matrix interfacial bond, as no mechancial clamping 
force (residual stress) is present. 
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confirming that this material had a weak interface. Note that varying the value of R n did not affect the stress level at 
which nonlinearity began, but instead only affected the stiffness and hardening slope of the tensile curve after 
debonding occurred. 


Effects of Residual Stresses 

The transverse tensile results presented to this point have included the effects of residual stresses. To further 
examine the effects of residual stresses on the transverse tensile results, analyses were conducted in which residual 
stresses were not applied. A four cell rectangular unit cell model with an R-ratio of 1.1 1 was used for the analyses, 
with the debond stress and interfacial compliance set to the values determined above. Simulated and experimental 
(ref. 26) stress-strain curves computed at a temperature of 650 °C (elevated temperature) are shown in figure 11, and 
results computed at 23 °C (room temperature) are shown in figure 12. As can be seen in the figures, at elevated tem- 
perature the application of residual stresses had minimal effect (as expected) on the computed results (particularly 
for strains less than 0. 1 percent), which are consistent with the results obtained for longitudinal tensile loading. 
Alternatively, at room temperature, when residual stresses were not included the stress level at which debonding 
occurred was significantly under predicted as compared to the experimental results. Also, the predicted stress values 
at high macroscopic strain levels are still somewhat higher than the experimental stress values, although the discrep- 
ancy is less than was seen at elevated temperature. This trend is consistent with the fact that (1) matrix inelasticity is 
less significant at room temperature as compared with elevated temperatures and (2) the interfacial model utilized 
does not allow for proper load shedding from the fiber to the matrix once debonding occurs. This fact will be dis- 
cussed more fully in the following discussion section. Clearly, these results verify once again that residual stresses 
need to be incorporated into all lower temperature composite simulations in order to obtain accurate results. 


Effects of Fiber Architecture 

Once the interface model was correlated, the influence of fiber packing on the transverse tensile response (due 
to loading along the 3-axis of fig. 3) of the SCS-6/TIMETAL 2 IS composite was examined for both elevated and 
room temperature. In these analyses, the normal debond stress and interfacial compliance were assumed to be inde- 
pendent of fiber packing, and the correlated values determined in the previous section were used. Five packing 
arrangements, which are displayed graphically in figure 13, are examined: square, square diagonal and three rectan- 
gular fiber packing (i.e., R = 0.74, R=l.ll (which is the average R ratio determined from micrographs for the 
actual material in (ref. 25) and not shown in fig. 13), and R = 1.34). The /?- ratio was defined as the ratio of the aver- 
age distance between fiber centers in the 3-direction (fig. 3) to the average distance between fiber centers in the 
2-direction. The square packing arrangement was equivalent to a rectangular pack with an tf-ratio of 1 .0. These 
fiber-packing architectures were chosen in order to be consistent with the qualitative analyses conducted in 
reference 6. 

Computed transverse tensile stress strain curves including residual stresses obtained at elevated temperature 
(650 °C) are shown in figure 14, and results computed at room temperature (23 °C) are shown in figure 15. In both 
cases, the computed results were compared to experimental results obtained from reference 26. At elevated tempera- 
ture, the fiber architecture did not appear to have a significant effect on the macroscopic stress level at which 
debonding occurred (indicated by the onset of nonlinearity ). However, the slope of the hardening curve was signifi- 
cantly affected by the fiber packing, with higher R-ratios yielding a softer response. The square diagonal pack 
resulted in the softest response, which was consistent with the qualitative results obtained by Arnold et al., (ref. 6). 
However, the very compliant post-yield slope is most likely an artifact of the unit cell model, since the square diago- 
nal unit cell has two fibers which can debond, thus leaving an intact cross-sectional area much smaller than the rect- 
angular unit cell model composed of a single fiber. Consequently, the remaining matrix will carry a greater stress 
and thus experience an increase in inelastic flow, thus producing a “softer” overall response history. To allow each 
unit cell model to have equivalent cross-sectional areas, a progressive debonding technique is required. While this 
capability was not present in MAC at the time these analyses were conducted, the concept has since been incorpo- 
rated into MAC/GMC and will be discussed in a future report. 

The results obtained at room temperature (see fig. 15) indicate once again that increasing the R-ratio results in a 
softening of the tensile stress-strain curve. Furthermore, once again the square diagonal pack yielded the softest 
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response. However, in contrast to the elevated temperature results, varying the fiber packing architecture at room 
temperature also affected the macroscopic stress level at which debonding occurred. These results illustrate that 
architecture influences residual stress build-up, so that as the R-ratio is increased, the stress at which debonding 
occurs is decreased. When the results at both temperature levels are examined simultaneously, the analytical results 
obtained by using a rectangular pack with an R ratio of 1 . 1 1 appears to provide a very good match to the experimen- 
tal values, which is reasonable since the average R ratio measured from micrographs is equal to 1.11. Overall, the 
results demonstrate that fiber packing has a significant effect on the overall response (and initial residual stress 
state), and the packing must be correctly modeled in order to correctly predict the material behavior. 

Arnold et al. (ref. 6) showed in a qualitative sense that using a more refined unit cell, which more completely 
approximated a circular fiber, yielded different results in simulations of transverse tensile tests than results com- 
puted using a square fiber approximation. Numerous other researchers have also examined the effects of fiber shape 
on the tensile response of composites, as reviewed and discussed in reference 6. For this study, the effects of varying 
the fiber shape, while maintaining the same fiber volume fraction and fiber packing were examined quantitatively 
through comparisons to experimental results. An R-ratio of 1.1 1 was used, and the interfacial compliance and 
debond stress values determined earlier were once again employed. A schematic of the unit cell model for the circu- 
lar fiber is shown in figure 16. An important point to note is that for the refined unit cell model, up to ten subcell 
interfaces may now debond once the debond stress reached. Loading was again applied along the 3-direction axis as 
defined in figure 3, and residual stresses were incorporated. 

Tensile curves obtained at 650° and 23 °C for both square and circular fiber shapes are shown in figures 1 7 and 
18, along with experimental results (ref. 26), respectively. The elevated temperature results indicated that varying 
the fiber shape primarily affected the slope of the hardening curve after debonding took place. However, the com- 
parison between the experimental and computed results was much better for the circular fiber model. This result was 
reasonable since the circular fiber model was more refined than the original square fiber model. The results com- 
puted using the circular fiber model at room temperature provided a better comparison to the experimental values in 
several respects. First of all, the original elastic modulus was predicted more accurately. Secondly, the slight break 
in the experimental curv e at a stress level of ~ 100 MPa was captured by the circular fiber model. The slight break in 
the experimental curve indicated that some debonding might have taken place at this point. Since the circular fiber 
unit cell allowed for the possibility of a gradual debonding, the model could capture this behavior. Furthermore, 
using the circular fiber model, which provided a better comparison to the experimental results, reduced the stiffness 
of the hardening curve. Overall, the results indicated that for modeling the transverse tensile response of a composite 
with a weak fiber/matrix interface, correctly modeling the fiber shape was very important. For the current study, 
however, one disadvantage of using the circular fiber shape was that the execution time for analyses using circular 
fibers increased by 20 to 30 times when compared to analyses conducted using square fibers. However, as discussed 
in reference 1 1, since the time that the analyses included in this report were conducted, significant modifications 
have been made to the MAC code which permits more computationally efficient analyses of complex fiber 
architectures. 


Discussion of Transverse Tensile Results 

For all of the computed transverse tensile results which have been presented, the three distinct regions of the 
stress-strain curve described by Majumdar and Newaz (ref. 3) were observed and predicted. Recall that these regions 
were denoted as follows. First, an initial linear region was present, wherein the macroscopic behavior remained elas- 
tic and the fiber and matrix remained bonded. A second linear region occurred next, wherein interfacial debonding 
occurred and the interface opened. Finally, a third nonlinear region occurred, in which significant matrix inelasticity 
was present. Furthermore, once the interfacial parameters were appropriately correlated and residual stresses prop- 
erly accounted for, the computed results compared reasonably well to the experimental values (ref. 26), particularly 
when a more refined unit cell model was utilized. However, the predictions made at elevated temperature generally 
over predicted the stresses in the inelastic portion of the tensile response, for reasons to be discussed below. 

To further analyze the reasons behind the over prediction of the transverse tensile stresses at elevated tempera- 
ture, the stress-strain response for each of the individual subcells was examined. In figures 19 and 20 the local 
stresses in the individual subcells (see fig. 3 for subcell numbering) are plotted; wherein subcell 1 corresponds to the 
fiber material and subcells 2, 3 and 4 are associated with the matrix material. In figure 19 the microstress verses 
strain response are shown with the residual strains present at the start of the mechanical loading (i.e., the curves do 
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not start at a strain value of zero, but at the actual residual strain value present after the cool-down process); 
whereas, in figure 20, the stress verses the total macroscopic strain in the composite is illustrated. In figure 20, for 
comparison the predicted total macroscopic transverse tensile stresses are also included in the plot. Furthermore, an 
analysis was conducted using MAC to simulate the tensile response of the bulk matrix. This analysis was conducted 
by assuming a negligibly small fiber volume fraction. The results from this additional analysis are also included in 
the figure (labeled as TIMET AL 2 IS). Clearly, all four subcells start off with a linear elastic response. Then at a 
microstress level of 96.5 MPa (the defined “debond stress”), the interface debonds, and the fiber subcell (subcell 1) 
and the adjoining matrix subcell (subcell 2) remain at a constant stress level. This constant stress is due to the nature 
of the interfacial debonding model utilized and the traction continuity that is imposed between subcells. The stress 
level in the two remaining matrix subcells (subcells 3 and 4) continue to increase, and in fact rise sufficiently to 
induce inelastic behavior in these subcells — thus explaining the onset of the third nonlinear region of the stress- 
strain curve as discussed above. However, realistically, one would expect the stress in the fiber and the adjoining 
matrix to reduce to zero subsequent to debonding; with the corresponding debond stress being transferred to the 
remaining matrix subcells after some time. This additional stress, that is transferred to the other matrix subcells, 
would most likely promote greater matrix inelasticity and consequently result in a further softening of the predicted 
macroscopic stress-strain curve; particularly at high temperatures where matrix inelasticity is quite significant. This 
deficiency in the interfacial model implemented within MAC is thought to be the cause of the over prediction of the 
macroscopic stress state in the composite. Note, one might question how is it that subcells 3 and 4 (which in this 
case are associated with the matrix material) can carry a microstress state higher than that of the matrix material 
itself as shown in figure 20. This is due to the fact that the subcell responses 3 and 4 depicted, are only one compo- 
nent of the multiaxial stress tensor. Consequently, if the effective (J 2 ) stress response of the matrix verses the macro 
strain response was plotted (not shown in the fig. 20) both should follow that of the matrix only response quite 
nicely. 

To further illustrate the multiaxial nature of the predicted stress-strain response of the [90] laminate at 650 °C, 
the variation of the macroscopic stress in the loading direction is plotted as a function of the macroscopic strain in 
each of the three coordinate directions in figure 21; that is, 1-1 (parallel to the fiber), 2-2 (perpendicular to the fiber 
and the loading direction), and 3-3 (perpendicular to the fiber and in the loading direction). Note that once again the 
symbol a is used to indicate graphically the approximate point at which the response deviates from proportionality. 
Clearly, t^e 1-1 response is essentially linear with only very small (Poisson’s ratio induced) strain levels being 
obtained. This response is reasonable and expected, as the response in this direction is dominated by the fiber, which 
is elastic and has a much higher stiffness than the matrix. Alternatively, in the transverse directions the responses 
shown indicate a linear region up to the on-set of debonding (Stage I), followed by another linear region at a reduced 
stiffness (Stage II), followed by significant inelasticity at high macroscopic stress levels (Stage III). 

A similar examination to that done above for the elevated temperature tensile test can be conducted for that of 
the response at room temperature. Again the micro stress-strain response (in the 33-loading, direction) of the indi- 
vidual subcells are plotted in figure 22 whereas, in figure 23 the corresponding subcell and composite stress states 
are shown versus the macro strain in the 33 direction. Once again, predicted results for the macroscopic composite 
and the TIMETAL 21S matrix are included for comparison. Similar to our previous elevated temperature observa- 
tions, once the interface debonds the fiber subcell (subcell 1) and its adjoining matrix subcell (subcell 2) remain at a 
constant stress level (the “debond stress”), since the present interfacial debond model incorporated does not allow' 
for any load shedding to the remaining matrix subcells. Here, however, the influence (development of residual stress 
and strain) of the initial cool-down is quite significant and apparent, in that (see figs. 22 and 23) there are significant 
residual strains in all four subcells (with the fiber subcell having the least amount of residual strain). Note the sig- 
nificant compressive residual stresses present within the fiber subcell and its adjoining matrix subcell (which 
explains the delay in interfacial debonding), and, the significant tensile residual stresses in the remaining two matrix 
subcells. Again, examining the variation of the macroscopic stress in the loading direction as a function of the 
macroscopic strain in each of the three coordinate directions, see figure 24, we see similar behavior to that observed 
at 650 °C but now the three stages of deformation are more distinct due to the presence of residual stresses. 


CROSS-PLY LAMINATE TENSILE RESPONSE 

To further validate the developed methodology, the tensile response of a [0/90] cross-ply laminate was simu- 
lated. A schematic of the unit cell employed is shown in figure 25. The fiber volume fraction, interfacial parameters 
and fiber architecture were the same as applied previously for the analyses with transverse tensile loading. Analyses 
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were again conducted at both 650 °C (elevated temperature) and 23 °C (room temperature). Residual stresses were 
applied through a cool-down process, and strain controlled tensile loads were subsequently applied along the 
3-direction axis (see fig. 25) at a strain rate of 1.0x1 0^/sec, in order to match experimental results obtained by 
Castelli (ref. 31). 

Simulated and experimental results (ref. 31) obtained at elevated temperature are shown in figure 26, and results 
obtained at room temperature are shown in figure 27. To demonstrate the effects of residual stresses, computed 
results in which residual stresses were not included are also shown in the figures. As shown previously, the incorpo- 
ration of residual stresses had minimal effects at elevated temperature on the resulting stress-strain response. At 
room temperature, however, the onset of nonlinearity was significantly effected by residual stresses. When residual 
stress effects were ignored, the onset of nonlinearity occurred at a much higher strain level (resulting in an over pre- 
diction of the stresses) as compared to the experimental results. When residual stresses were incorporated into the 
analysis, a much better match with the experimental results was obtained. At elevated temperature, there was some 
discrepancy between the experimental and computed results, which could have been due to the interfacial model 
formulation, which did not permit load shedding once debonding occurred as discussed earlier in this report. Fur- 
thermore, fiber damage prior to final failure might have occurred in the [0°] plies in the experimental tests which 
was not accounted for in the model (as discussed earlier). 


CONCULSIONS 

The effects of residual stresses, fiber packing and shape, and the influence of interfacial bonding on the tensile 
stress strain response of a titanium matrix composite (SCS-6/TIMETAL 2 IS) has been modeled for a composite 
under uniaxial longitudinal and transverse loads, as well as for a cross-ply laminate. Additionally, the analytical 
results have been compared to experimentally obtained results in order to provide a quantitative characterization of 
the various parameters examined. For all of the laminates studied, proper application of residual stresses was 
required in order to accurately model the composite response. Furthermore, the fiber architecture needed to be prop- 
erly described and the fiber/matrix interface needed to be appropriately characterized in order for the analytical 
response to compare favorably to the experimental results. Through proper characterization of the fiber/matrix inter- 
face model described in this work, the stress level at which debonding occurred and the slope of the hardening curve 
were appropriately simulated. These analyses thus indicated that the generalized method of cells micromechanics 
model in combination with the generalized viscoplastic potential structure inelastic constitutive model could 
adequately predict the actual behavior of the material. 

Future work will encompass several points. First, as discussed earlier, the interfacial model utilized in this 
report, while easy to utilize and computationally efficient, has significant deficiencies, particularly in terms of prop- 
erly shedding load from the fiber to the matrix once debonding has occurred. A need is present to develop an inter- 
face model where the interface is modeled as a distinct constituent, with material properties that are different in 
tension and compression (to account for application of residual stresses), and allows for the proper stress transfer 
between the fiber and surrounding matrix. Additionally, creep and relaxation analyses, as well as cyclic analyses, 
should be conducted in order to examine the effects of the parameters studied in this report for loadings other than 
monotonic tensile loadings. 
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TABLE 1.-TT-2IS MATERIAL PARAMETERS 



23 °C 

300 °C 





1 14. 1 


95.1 


59.7 

CTE (p 0 C"' ) 

7.717 

9.209 

10.70 

12.13 

14.09 

k ( MPa) 

1029 

768.4 

254.2 

5.861 

0.756 

p (MPa/sec) 

667.6 

137.8 

1.45x10'’ 

6.19x10^ 

1.13x10'" 

£ 0 (MPa) 

6.908 xlO' 5 

1.035x1 O’* 

2.756x10'' 

5.870x1 O' 1 

6.346x10'' 

R a ( 1/sec) 

0 

0 

1.68x1 O' 7 

l.OOxlO' 6 

6.00x1 0" 5 

p 

0.001 

0 

0 

0 

0 


Temperature-independent: v = 0,365, n = 3.3, B l - 0.05, p- 1.8, q - 1.35 


Double 

periodicity 


<1 <2 




P = N p 


P = 2 
P = 1 


y=1 7=2 


h N P 


h 2 

h 1 t 


7 = N 7 



Figure 1 . — Generalized Method of Cells, subcell dimension nomenclature for both (a) double and (b) triple 
periodicity conditions (reference 8). 



Figure 2. — Micrograph of SCS-6/TIMETAL 21 S composite specimen. 
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Figure 3. — Two-dimensional GMC model utilized to 
simulate longitudinal and transverse tensile loading 
(reference 8). For longitudinal loading, the loading dir- 
ection was along the 1 direction axis. For transverse 
loading, the loading was along the 3 direction axis. 



Strain, % 


Figure 4.— Predicted longitudinal tensile stress-strain response of SCS-6/TIMETAL 21 S at 650 °C 
with comparison to experiment (reference 27). 
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Figure 5. — Predicted longitudinal tensile stress-strain response of SCS-6/TIMETAL 21 S at 23 °C 
with comparison to experiment (reference 27). 



°0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 


Strain, % 

Figure 6. — Predicted longitudinal tensile stress-strain response of 
SCS-6/TIMETAL 21 S at 23 °C with fiber response and fiber failure 
stresses Included. 


NASA/TM— 2000-209798 


16 





Stress, MPa 




°' / / 

t// 

v 


O Experimental 

ctqq = 0 MPa 

cr DB = 34.48 MPa 

— - — <TQg = 68.95 MPa 
°bB = ®6.53 MPa 
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Figure 7. — Effect of debond stress on the predicted transverse tensile response of 
SCS-6/TI METAL 21 S at 650 °C with comparison to experiment (reference 26). 
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°DB = 34 48 MPa 
o-qb = MPa 
<r n o = 96.53 MPa 


3.0 0.2 0.4 0.6 0.8 1.0 1.2 

Strain, % 

Figure 8. — Effect of debond stress on the predicted transverse tensile response of 
SCS-6/TIMETAL 21 S at 23 °C with comparison to experiment (reference 26). 


NASA/TM — 2000-209798 








Stress, MPa 
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Figure 9. — Effect of interfacial compliance on the predicted transverse tensile response of 
SCS-6/TIMETAL 21 S at 650 °C with comparison to experiment (reference 26). 



Figure 10. — Effect of interfaciai compliance on the predicted transverse tensile response of 
SCS-6/TIMETAL 21 S at 23 °C with comparison to experiment (reference 26). 
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Strain, % 

Figure 1 1 . — Effect of residual stresses on the predicted transverse tensile response of 
SCS-6/TIMETAL 21 S at 650 °C with comparison to experiment (reference 26). 
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Figure 12. — Effect of residual stresses on the predicted transverse tensile response of 
SCS-6/TIMETAL 21 S at 23 °C with comparison to experiment (reference 26). 
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Figure 1 3. — Schematics of representative volume elements of the 
fiber packing arrays examined, including square, rectangular with 
R = 0.74 and R = 1 .34, and square diagonal (reference 8). 
Transverse tensile loading is along 3 direction axis. 
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MPa -n Stress, MPa 



igure 14. — Effect of fiber packing arrangement on the predicted transverse tensile response of 
SCS-6/TIMETAL 21 S at 650 °C with comparison to experiment (reference 26). 
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Figure 15. — Effect of fiber packing arrangement on the predicted transverse tensile response of 
SCS-6/TIMETAL 21 S at 23 °C with comparison to experiment (reference 26). 
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Circular fiber 


Figure 16. — Schematic of representative volume 
element with circular fiber shape (reference 8). 
Transverse tensile loading is along 3 direction 
axis. 
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Figure 17. — Effect of fiber shape on the predicted transverse tensile response of SCS-6/TIMETAL 
21 S at 650 °C with comparison to experiment (reference 26). 
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Macro strain, % 


Figure 20. — Predicted transverse tensile stress vs. macroscopic strain response of individual 
subcells of simulated SCS-6/TIMETAL 21 S composite at 650 °C. Predicted macroscopic re- 
sponse of composite and predicted response of bulk TIMETAL 21 S included for comparison. 



Micro strain, % 


Figure 21. — Predicted transverse tensile stress/strain response of simulated SCS-6/TIMETAL 
21 S composite at 23 °C with off-axis strains included. 
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Figure 22. — Predicted transverse tensile stress vs. micro strain response of individual subcells of 
simulated SCS-e/TIMETAL 21 S composite at 23 °C. 
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Figure 23. — Predicted transverse tensile stress vs. macroscopic strain response of individual 
subcells of simulated SCS-6>TIMETAL 21 S composite at 650 °C. Predicted macroscopic re- 
sponse of composite and predicted response of bulk TIMETAL 21 S included for comparison. 
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Strain, % 


Figure 24. — Predicted transverse tensile stress/strain response of simulated SCS-6/TI METAL 
21 S composite at 23 °C with off-axis strains included. 



Figure 25. — Three-dimensional GMC model used to simulate tensile 
response of SCS-6/TIMETAL 21 S [0/90] laminate (reference 8). 
The loading direction was along the 3 direction axis. 
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Figure 26. — Predicted tensile stress-strain response of SCS-6/TIMETAL 21 S [0/90] laminate at 
650 °C with comparison to experiment (reference 31). 
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Figure 27. — Predicted tensile stress-strain response of SCS-6/TIMETAL 21 S [0/90] laminate at 
23 °C with comparison to experiment (reference 31). 
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